Simulated nuclear spin-lattice relaxation in Heisenberg ferrimagnets: 
Indirect observation of quadratic dispersion relations 



Shoji Yamamoto 

Department of Physics, Okayama University, Tsushima, Okayama 700-8530, Japan 

(26 October 1999) 

In response to recent proton spin relaxation-time measurements on NiCu(pba)(H20)3-2H20 with 
pba = 1, 3-propylenebis(oxamato), which is an excellent one-dimensional ferrimagnetic Heisenberg 
model system of spin-(l, |), we study the Raman relaxation process in spin-(S', s) quantum ferri- 
magnets on the assumption of predominantly dipolar hyperfine interactions between protons and 
magnetic ions. The relaxation time Ti is formulated within the spin-wave theory and is estimated 
as a function of temperature and an applied field If by a quantum Monte Carlo method. The 
low-temperature behavior of the relaxation rate T^ 1 qualitatively varies with (S,s), while T^ 1 is 
almost proportional to H -1 ^ 2 due to the characteristic dispersion relations. 

PACS numbers: 75.10.Jm, 75.40. Mg, 75.30. Ds, 76.60.Es 



Quantum spin chains have been providing broadband 
topics and have indeed invited numerous researchers 
to explore into them. One of the hot topics is the 
mixed-spin system with two kinds of antiferromagnet- 
ically exchange-coupled centers. Besides the competi- 
tion between the massive and massless states Q which 
is an interesting issue in itself, quantum ferrimagnetic 
phenomena have been attracting wide interest. There 
already exists an accumulated chemical knowledge on 
bimetallic chain compounds. Kahn et at j^j succeeded 
in synthesizing the variety of ferrimagnetic bimetallic 
chains such as MCu(pba)(H 2 0) 3 -nH 2 (M = Mn,Ni) 
g and MCu(pbaOH)(H 2 0) 3 -nH 2 (M = Fe,Co,Ni) §), 
where pba — 1, 3-propylenebis(oxamato) and pbaOH = 
2-hydroxy-l,3-propylenebis(oxamato). Under close con- 
tact between chemists and physicists, not only structural 
but also magnetic properties of them were elucidated || . 
It may be a few mathematical investigations [0-0] that 
renewed our enthusiasm toward the subject. Motivated 
by these works, several authors carried out spin- wave 
analyses [p 10 on the alternating-spin Heisenberg chains 
and pointed out the coexistence of the ferromagnetic and 
antiferromagnetic branches in the low- lying excitations. 
Modified spin- wave theories Jo|l2|] were further applied 
to the thermodynamics [|13|]14[ | and a ferromagnetic- to- 
antiferromagnetic crossover was found in the tempera- 
ture dependences of the specific heat and the magnetic 
susceptibility. The dual structure of the excitations re- 
sults in rich physics in a magnetic field as well, e.g. the 
double-peak structure of the specific heat jljj and quan- 
tized plateaux in the ground-state magnetization curve 
|fL6| , [L7f . Quantum ferrimagnetism and related phenom- 
ena are still central subjects in the field of material sci- 
ence. A family of one-dimensional oxides, SrsMM'Og 
(M = Ni, Cu, Zn; M' = Pt,Ir), consisting of alternat- 
ing M'C>6 octahedra and MOq trigonal prisms, visualizes 
the competition between various magnetic phases fll8| , 
whereas an ordered double perovskite, Sr 2 FeMoC>6, acts 
as a half-metallic ferrimagnet |l9| , ^o| ; exhibiting a novel 
magnetoresistive behavior. 



In such circumstances, the proton spin relax- 
ation time T\ has recently been measured (^] for 
NiCu(pba)(H 2 0)3-2H 2 0, which is typical of the spin- 
(1, i) ferrimagnetic Heisenberg chain, in an attempt to 
reveal the electron-spin dynamics inherent in quantum 
ferrimagnets. Although the susceptibility-temperature 
product, %T, shows a minimum in its temperature depen- 
dence (see Fig. || bellow), the observed relaxation rate 
Tj -1 is monotonically decreasing with the increase of tem- 
perature. The most interesting observation is the depen- 
dence of the relaxation rate on the applied magnetic field 
H, which looks like oc if -1 / 2 , though the authors did 
not rule out a logarithmic behavior, T^ 1 ~ lnif, based 
on a rather approximate argument. We here present 
a rapid communication on our numerical simulation of 
T\ all the more because the measurements reported are 
highly stimulative but still preliminary. We show that T\ 
should indeed act as T^ 1 oc H^ 1 / 2 for general spin-(S, s) 
ferrimagnetic Heisenberg chains due to the characteristic 
dispersion relations. The if ~ ^-dependence of T^ 1 may 
first remind us of the diffusive behavior of the spin cor- 
relation function p2| . However, a totally different mech- 
anism may cause the field-dependent relaxation rate in 
the present system. Although the dynamic behavior |2^| 
of quantum ferrimagnets is a fascinating subject, mea- 
surements in this direction seem to be still in their early 
stage. We expect the present calculation to accelerate 
experimental investigations into the dynamic properties. 

We consider two kinds of spins S and s (S > s) alter- 
nating on a ring with antiferromagnetic exchange cou- 
pling between nearest neighbors, as described by the 
Hamiltonian 

N N 

H = JJ2 (Sj ■ Sj + Sj ■ S j+1 ) - g^H^S^ + s z 3 ) , 

(1) 

where we have set the g factors of spins S and s both 
equal to g. Keeping in mind ferrimagnetic bimetallic 
chain compounds such as NiCu(pba)(H 2 0)3-2H 2 M, we 
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FIG. 3. Field dependences of the proton spin relaxation 
rate at various values of temperature: (a) (S,s) — (1, ~), (b) 
= (§,§), and (c) (S, a) = (§,!). 
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+ (A z sinh 2 e k - a z cosh 2 6 k fn+{n+ + 1)] , (6) 



where A z = Y^j A z and a z = J^j a j are the = 
Fourier components of the hyperfme coupling constants, 
and nj~ = (ata^) and = {0\,(3 k ) are the thermal aver- 
ages of the numbers of the spin waves at a given temper- 
ature. We then have an idea of evaluating n k by a quan- 
tum Monte Carlo method. Though the divergence of the 
ground-state magnetization, which plagues the spin- wave 
treatment of low-dimensional antiferromagnets, does not 
persist in our model, n^r still diverge as temperature in- 
creases as far as we naively estimate them with the non- 
compact bosonic Hamiltonian (||). However, introducing 



a certain constraint, for instance, on the magnetization, 
we can totally remove the difficulty and obtain a fine de- 
scription of the thermodynamics jl3| , where it should be 
stressed that the dispersion (|J) is still preserved. Thus, 
relying upon the Bogoliubov transformation, we perform 
the Monte Carlo sampling for the relevant spin operators 



= coahfc^ fc(j ._ i)sf ^m,y eikU+i 
V2N ^ 3 y/N t-r 1 



).«+ 



4^ 3 Vn ^ 



(7) 



with the original compact Hamiltonian ([!]). While we 
base the relaxation process on the spin- wave excitations, 
we take grand-canonical averages within the original sys- 
tem. The thus-estimated are indeed consistent with 
the exact calculations [23] of the dynamic structure fac- 

tors S ± (k,cu) = £„|<n|Sfc +*fc|g)f - E n + E s ) 
through the relation = 5"^ (fc, ^) , where |g) is the 
ground state with energy E g and magnetization (S—s)N. 
Since the applied field H is small enough in practice 
(gH^H < 10 -2 J) | pl[ , we neglect the Zeeman term of the 
Hamiltonian (jl|) in the numerical treatment. Considering 
the case of the contributing protons being located near 
the smaller-spin magnetic ions ||, we set A z ja z equal to 
1/3. 

We show in Fig. [I] the relaxation rate as a func- 
tion of temperature. The calculations in the case of 
(S,s) = (1, i) are qualitatively consistent with the ex- 
perimental observations |2l| ], where T^ 1 is a monoton- 
ically decreasing function of temperature. T^ 1 hardly 
depends on temperature except for the effectively-low- 
temperature region satisfying k^T < W~ , where W~ is 
the ferromagnetic band width, which may qualitatively 
be identified with cj r =7r — ^ r =0 . Considering the dis- 
persion relations (U), it is convincing that the ferromag- 
netic, rather than antiferromagnetic, nature dominates 
the low-temperature relaxation process. In this context, 
we show in Fig. |^ the numerical calculations of \T , 
which is closely related with Ty 1 |2f| . %T is a monotoni- 
cally decreasing function in ferromagnets, while a mono- 
tonically increasing function in antiferromagnets. Thus 
we learn that ferrimagnets display the mixed behavior 
in their temperature dependences of yT . Recently it 
has been pointed out [P6| that Spin-(5, s) ferrimagnetic 
chains behave similar to combinations of spin-(S' — s) fer- 
romagnetic and spin- (2s) antiferromagnetic chains pro- 
vided S — 2s, which is again illustrated in Fig. |^. For 
S > 2s, quantum ferrimagnets are predominantly ferro- 
magnetic rather than antiferromagnetic, where \F more 
rapidly decreases than the, let us say, balanced behav- 
ior, ( x (S-s)-fcrro + x (2 S )-antifcrro- )Tj w h erea s for S < 2s, 

vice versa, where the low-tempcrature decrease is duller 
than the balanced behavior. It is quite interesting to ob- 
serve Fig. [j] from this point of view. The relatively rapid 
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decrease of Tf 1 at low temperatures for (S,s) = (|, |) 
must be of predominantly ferromagnetic aspect, while 
the shoulder-like, even increasing, behavior at low tem- 
peratures for (S,s) = (§,1) can be regarded as an anti- 
ferromagnetic feature. However, the overall gradual in- 
crease at high temperatures, as observed for \T , does 
not appear. It is the field-dependent prefactor to nj: in 
Eq. (jfj), coming from the energy-conservation require- 
ment 5(E m — E n -wn) in Eq. (||), that suppresses the 
latent increasing behavior at high temperatures. Since 
cjjm <C J, the k = component is predominant in the 
summation in Eq. (^). In general, with the increase 
of temperature, decrease in the vicinity of the zone 
center, whereas they increase near the zone boundary 
p3| , where we note that both have their peaks at 
k = 0. Therefore, the increase of the field, reducing the 
predominance of the k = component, smears out the 
temperature dependence of T-f . 

The characteristic prefactor in Eq. (^J) further leads to 
the unique field dependence of T^ 1 . In Fig. || we plot 
T^ 1 as a function of = TN-ff- T± is highly linear 
with respect to if -1 / 2 . We stress that the observations in 
Fig. U originate from the quadratic dispersion relations 
peculiar to this system and therefore the H — > oo ex- 
trapolation of T x ~ results in a considerably small value. 
Thus the present phenomena are distinguishable from the 
diffusive behavior, (St(t)SZ„(0)) oc c~ Dq which also 



causes the field dependence of the form |22[ 



^ = P + Q\ J jf 



(8) 



The constant P is inherent in the fluctuation-dominant 
spin dynamics. For (CHs^NMnCkj and L1V2O5, which 
are both well known to be Heisenberg linear-chain com- 
pounds whose spin correlation functions behave diffu- 
sively at long times, P/Q was estimated to be 139 
p2f (in the high-temperature limit) and 376 |^] (at 
k-^T/J ~ 0.65), respectively. On the other hand, 
the spin-(l, i) Heisenberg ferrimagnetic chain compound 
NiCu(pba)(H 2 0) 3 -2H 2 exhibits P/Q ~ 12 [|lj (at 
k-QT/J ~ 2.3), which is much smaller than those of the 
spin-diffusive materials but is in good agreement with 
the present numerical findings (Table ||) . Therefore the 
T\ measurements could implicitly demonstrate the low- 
energy quadratic dispersion relations of quantum ferri- 
magnets. When the ATY-type exchange anisotropy moves 
the model away from the Heisenberg point, the chain 
turns critical, showing a linear dispersion |rj|]28||, and the 
present field dependence of T^ 1 should no more be ex- 
pected. Thus the nuclear-magnetic-resonance, as well 
as neutron-scattering, measurements allow us to reveal 
the low-energy structure efficiently. We believe that the 
present calculations will greatly motivate further experi- 
mental investigations into quantum ferrimagnets. 
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comments and fruitful discussion on proton-NMR mea- 
surements for quantum ferrimagnets. This work was sup- 
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TABLE I. P/Q in the case of Ti being fitted to 
Tf 1 = P + Q{^H/jy 1/2 . The estimates in the limit of 
UbT / J — ► oo are obtained through a rough extrapolation and 
are thus no more than for reference. 
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